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14.  ABSTRACT  - - - 

The  over-arching  objective  of  this  research  was  to  improve  the  quantitative  basis  for  prediction  of  the  occurrence  of  Coronal  Mass 
Ejections  (CMEs)  and  their  impact  on  space  weather.  The  most  directly  useful  discovery,  from  the  point  of  view  of  Space  Weather 
forecasting,  is  that  a  CME  associated  with  an  erupting  X-ray  sigmoid  is  more  likely  to  have  flux  rope  (“magnetic  cloud”)  structure 
at  1  AU  than  an  arbitrarily-selected  CME.  On  impacting  the  magnetosphere,  the  sigmoid-associated  flux  ropes  in  our  database  all 
produced  at  least  "moderate"  (Dst  >  50  nT)  geomagnetic  storms,  and  several  produced  major  storms.  Hence,  not  only  are  sigmoidal 
active  regions  more  likely  to  erupt,  as  we  showed  in  earlier  work  (Canfield  Hudson  &  McKenzie  1 999),  they  are  more  likely  to 
cause  significant  storms  when  they  do.  Also  discovered  was  a  signature  of  impending  eruption  in  H-alpha  imaging  spectra  of 
active  regions,  which  could  be  used  to  predict  the  occurrence  of  eruptions  a  few  hours  in  advance,  whereas  the  sigmoid  signature  is 
more  useful  on  time  scales  of  days.  These  discoveries  have  been  documented  in  1 1  publications  and  26  presentations. 
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OBJECTIVES 

•  Improve  the  quantitative  basis  for  prediction  of  the  occurrence  of  Coronal  Mass  Ejections  (CMEs) 
and  their  impact  on  space  weather. 

•  Extend  earlier  work  into  the  regime  of  quantitative  prediction  based  on  vector  magnetograms,  X-ray 
images,  and  other  data 

•  Learn,  through  3D  MHD  models,  why  active  regions  that  are  large,  or  show  twisted  fields,  are  the 
most  likely  sites  of  CMEs. 

•  Understand  the  form  and  characteristics  of  current  singularities  in  2D  and  3D  steady  states,  with 
emphasis  on  the  effects  of  helicity. 

ACCOMPLISHMENTS/NEW  FINDINGS 

Data  Analysis 

We  developed  a  database  of  eruptive  solar  events  that  includes  Mees  Solar  Observatory  active  region  vector 
magnetograms,  Yohkoh  X-ray  sigmoids,  cusps  and  jets,  interplanetary  data  from  ACE,  WIND  and  IMP,  and 
geomagnetic  indices.  We  used  the  database  to  study  solar,  interplanetary,  and  geophysical  effects  of  erupting 
solar  active  regions.  With  this  database  we  made  new  findings: 

•  A  CME  associated  with  an  erupting  coronal  sigmoid  is  more  likely  to  have  flux  rope  (called 
“magnetic  cloud”)  structure  at  1  AU  than  an  arbitrarily-selected  CME  (Leamon,  Canfield,  and 
Pevtsov,  2002).  On  impacting  the  magnetosphere,  the  flux  ropes  in  our  database  all  produced  at  least 
"moderate"  (Dst  >  50  nT)  magnetic  storms.  This  result  has  importance  in  Space  Weather 
forecasting:  it  shows  that  not  only  are  sigmoidal  active  regions  more  likely  to  erupt,  as  we  showed  in 
earlier  work  (Canfield  Hudson  &  McKenzie  1999),  they  are  more  likely  to  cause  significant  storms 
when  they  do. 
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•  The  correlations  between  coronal  and  1  AU  properties  of  sigmoidal  eruptions  (most  significantly, 
leading  magnetic  field  direction)  reverse  at  solar  minimum,  when  sunspot  polarity  reverses,  as 
opposed  to  solar  maximum,  when  the  global  solar  dipole  reverses  (Pevtsov  &  Canfield,  2001).  In 
other  words,  eruptions  associated  with  active  regions  do  not  follow  the  general  rule,  derived  without 
regard  to  active  region  association,  which  says  that  the  global  solar  dipole’s  orientation  solely 
determines  the  leading  polarity  of  the  interplanetary  magnetic  cloud. 

•  The  overall  twist  of  active  region  magnetic  fields  is  very  small  compared  to  that  of  interplanetary 
magnetic  clouds  that  are  associated  with  them  (Leamon  et  al  2003,  in  preparation).  Roughly 
speaking,  the  helical  structure  of  magnetic  clouds  inferred  from  Bessell-function  fits  to  the 
interplanetary  data  contains  more  than  order  of  magnitude  more  turns  than  the  sheared  structure  of 
active  regions  inferred  from  vector  magnetogram.  This  discovery  is  important  because  the  only 
models  of  CMEs  that  can  explain  this  result  are  those  that  invoke  magnetic  reconnection  to  form  the 
magnetic  cloud  from  the  pre-eruption  corona. 

•  Yohkoh  SXT  image  sequences  of  the  shapes  of  191  X-ray  sigmoids  and  the  histories  of  eruption 
(evidenced  by  cusp  and  arcade  signatures)  of  their  associated  active  regions  reveal  that  the 
distribution  of  sigmoid  shapes  is  quite  narrow,  confirming  an  earlier  result  by  Rust  and  Kumar 
(1996),  and  the  frequency  of  eruption  does  not  depend  significantly  on  shape  (Leamon  et  al.,  2003). 
We  know  of  no  explanation  of  this  interesting  clue  about  the  origin  of  this  phenomenon. 

•  Vector  magnetograms  give  values  of  the  large-scale  total  twist  of  active  regions  in  which  flare- 
related  Yohkoh  SXT  signatures  of  eruption  are  observed  (Leamon  et  al,  2003).  We  find  no  evidence 
of  eruption  for  values  of  total  twist  remotely  approaching  the  threshold  for  the  kink  instability.  In 
other  words,  we  find  no  evidence  that  the  MHD  kink  instability  plays  a  significant  role  in  solar 
eruptions. 

•  We  discovered  a  signature  of  impending  eruption  in  H-alpha  imaging  spectra  of  active  regions.  This 
signature  is  present  during  a  three-hour  period  before  flares  that  are  eruptive,  but  is  absent  prior 
those  that  are  not  (Des  Jardins  and  Canfield,  2003).  This  is  a  statistically  significant  result,  based  on 
a  study  of  16  events.  These  properties  could  be  used  to  predict  the  occurrence  of  eruptions  a  few 
hours  in  advance,  whereas  the  sigmoid  signature  is  more  useful  on  time  scales  of  days. 

Finally,  we  wrote  a  review  of  sigmoids  and  eruptions  (Canfield,  Hudson,  and  Pevtsov  2000)  for  non-solar 
space  weather  researchers. 

Theory  and  Modeling 

We  completed  a  3-dimensional  numerical  MHD  simulation  of  the  emergence  of  twisted  flux  tubes  into  the 
corona.  The  simulation  was  initialized  with  a  straight  tube  of  twisted  magnetic  field  located  in  the  upper 
convection  zone.  Buoyancy  drives  an  arched  segment  of  a  twisted  flux  tube  upward  through  the 
photospheric  layer  and  into  the  corona.  Matter  drains  from  the  coronal  field,  which  thereafter  undergoes  a 
rapid  expansion  and  attains  sigmoidal  topology  (Maraga  and  Longcope  2001).  We  discovered  that  magnetic 
field  lines  that  emerge  with  a  large  aspect  ratio  (the  ratio  of  height  to  footpoint  separation)  simply  continue 
to  expand  into  the  outer  atmosphere,  while  field  lines  emerging  with  a  small  aspect  ratio  do  not  (Magara 
2001).  This  is  a  previously  unknown  aspect  of  the  evolution  of  active  regions. 

We  analyzed  the  injection  of  magnetic  energy  and  magnetic  helicity  into  the  atmosphere  during  the 
simulated  flux  emergence  (Magara  &  Longcope  2003).  We  found  that  the  emergence  contributions  are  the 
dominant  ones  at  the  early  phase  of  flux  emergence  and  that  later  that  role  is  played  by  the  shearing 


contributions.  The  emergence  starts  with  a  simple  dipole  structure  formed  in  the  photosphere,  which  is 
subsequently  deformed  and  fragmented,  leading  to  quadrupolar  structure. 


We  derived  a  new  mathematical  scheme  for  mapping  the  connectivity  of  a  potential  magnetic  field  arising 
from  an  arbitrary  distribution  of  discrete  sources  (Longcope  &  Klapper  2002).  The  field  lines 
interconnecting^  sources  are  classified  into  domains,  defining  the  field's  connectivity.  The  sche™ 
enables  one  to  quantify  how  much  energy  is  stored  in  active  regions,  and  where  it  is  stored  in  topological 

terms. 
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